Draft version June 8, 2012 

Preprint typeset using I^'T^]X style emulateapj v. 5/2/11 



PRECISION COSMOGRAPHY WITH STACKED VOIDS 

GUILHEM LAVAUX"''^"'^''^ 

Department of Physics, University of Illinois at Urbana-Cliampaign, 1002 W Green St, Urbana, IL, 61801, USA 



Benjamin D. Wandelt"'-"'''' 

UPMC Univ Paris 06, UMR 7095, Institut d'Astrophysique de Paris, 98 bis, boulevard Arago, 75014 Paris, France 

Draft version June 8, 2012 

ABSTRACT 

We present a purely geometrical method for probing the expansion history of the Universe from the 
observation of the shape of stacked voids in spectroscopic redshift surveys. Our method is an Alcock- 
Paczyhski test based on the average sphericity of voids posited on the local isotropy of the Universe. 
It works by comparing the temporal extent of cosmic voids along the line of sight with their angular, 
spatial extent. We describe the algorithm that we use to detect and stack voids in redshift shells on 
the light cone and test it on mock light cones produced from iV-body simulations. We establish a 
robust statistical model for estimating the average stretching of voids in redshift space and quantify 
the contamination by peculiar velocities. Finally, assuming that the void statistics that we derive 
from A'^-body simulations is preserved when considering galaxy surveys, we assess the capability of 
this approach to constrain dark energy parameters. We report this assessment in terms of the figure 
of merit (FoM) of the dark energy task force and in particular of the proposed EUCLID mission which 
is particularly suited for this technique since it is a spectroscopic survey. The FoM due to stacked 
voids from the EUCLID wide survey may double that of all other dark energy probes derived from 
EUCLID data alone (combined with Planck priors). In particular, voids seem to outperform Baryon 
Acoustic Oscillations by an order of magnitude. This result is consistent with simple estimates based 
on mode-counting. The Alcock-Paczyhski test based on stacked voids may be a significant addition 
to the portfolio of major dark energy probes and its potentialities must be studied in detail. 



1. INTRODUCTION 



The physical nature of Dark Energy, detecte d through 
supcrnovae luminosity distance measurements ( [Kowalski 



et al. 2008), Baryonic Acoustic Oscillations (BAo, Per 
aval ct al. 2010| and t he Cosmic Microwave Background 
i jKora atsu ct ah 201 1|), st ill evades us. The BOSS sur- 
vey fechlcgcl et al. '5007) has been designed to assess 
whether the equation ot state of Dark Energy is indeed 
constant and equal to minus one. However, observations 
based on baryonic acoustic oscillations are limited by the 
minimal volume required to estimate the scale of these 
oscillations, typically '-^lOO/i"^ Mpc. With the advent 
of large galaxy spectrosco pic redshift survey, su ch as the 
Sloan Digital Sky Survey ( [Abazajian et al.|2009| ) , we now 
have access to a three-dimensional representation of the 
large-scale structure on our light-cone on vastly different 
scales. 

The we l l kno wn Alcock-Paczyhski test (Alcock & 



Paczynski||1979 l can be applied to any structure tor 
which we know the physical size or, more weakly, the ra- 
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tio of its extent along the line of sight and its angular size. 
In particular, if we had a population of standard spheres 
scattered throughout cosmic history we could measure 
the cosmological expansion directly. Absent such a pop- 
ulation, the next best thing is a population of objects 
whose average shape is spherical. 

Cosmic voids are such a population and hence promis- 
ing candidates for probing the expansion geometry of the 
Universe. Even though individual void shapes may be 
complicated, the average void is spherical in an isotropic 
and homogeneous universe. Detecting all voids observed 
in a galaxy survey and stacking voids of similar sizes 
and redshifts projects out the details of individual void 
shapes. Since the average shape of voids is known to be 
spherical, the observed, stretched shape in redshift space 
is a direct function of the local Hubble expansion of the 
Universe and the angular diameter distance at the red- 
shift of the void and hence a sensitive function of the 
cosmological parameters, in particular those parameter- 
izing the dark energy equation of state. 

Voids are spatially localized coherent structures with 
sizes between ~5-20 Megaparsecs. If voids can be used 
to construct an Alcock-Paczyhski test, they will unlock 
a much larger number of modes for precision cosmol- 
ogy than are accessible by the BAO technique, including 
modes that are in the mildly non-linear regime. 

Crucially, by focusing on void regions we take advan- 
tage of the much more easily modeled phase-space struc- 
ture in low-density regions compared to high-density re- 
gions. 

The promise of an AP test based on voids was first no- 
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ticed by Ryden ( 1995 1 who proposed using the apparent 
stretching o± void shapes in redshift space coordinates 
to estimate the local geometry of the expansion. Doing 
so properly requires selecting voids that have the same 
overall size and density. This work has then yielded a 



oscillations (Ballinger et al. 11996 


Ryden & Melott||1996 


Schmidt et ai. 2UU1), though the complicated shapes o 



logical information a t high signa l-to-noise. 
In a similar spirit, Jimenez & Locb (2002 



test which uses cosmic clocJts as tracers ot 



I proposed a 
cosmic time 



to which we can compare the measured galaxy redshifts. 
The Alcock-Paczyhski test can be thought of as a differ- 
ential version of this approach — one could say comparing 
the radial, temporal extent of voids with their angular, 
spatial extent amounts to using them as "cosmic stop- 
watches." Further, it does not rely on spectral modeling 
t o extract galaxy age s. 



Park & Kim ( 20To| proposed using genus statistics to 
measure expansion, 'i'heir claim was that the genus is in- 
sensitive to redshift space distortions while pe culiar ve- 
locities will mildly affect the void technique (Ballinger 
1996|. However, the genus is not a spatially lo- 



et al. 



cal quantity, which may be prone to a number of obser- 
vational problems like non-trivial edge effects and inho- 
mogeneous incompleteness corrections. It is also model 
dependent. 

Our method requires spectroscopic survey since red- 
shift errors in photometric redshift catalogs wash out the 
line-of-sight information on the scale of all but th e most 
extreme vo ids even with ~ 0.7% precision (e.g. Ilbert 
et al.|2009 ). An order of magnitude improvement in red- 
shift precision would likely be required to directly observe 
any non-linear th ree-dimensional structures (c/. [Jasche 
& Wandelt (2011) for a possible approach). 

Our paper is organized as follows. In Section [2j we 
show how voids may give us a direct probe of cosmolog- 
ical parameters through shape stretching. In Section |3j 
we explain the method we use for finding and stacking 
voids on the light-cone and inferring the local expansion 
from the shape of stacked voids. In Section|4| we test our 
method on TV-body simulations. We derive the profiles of 
stacked voids in simulation, the sensitivity to contamina- 
tion of the redshifts by peculiar velocities and the number 
density of voids found in the simulation. In Section [s] 
we derive the Hubble-diagram of expansion in the Simula^ 
tion from our mock-observation of voids. We do a Fisher- 
matrix analysis of the measurement of Dark Energy phys- 
ical parameters using the expansion rate derived from 
voids. We apply this formalism to the survey specifi- 
cations of the main ga laxy sample of the Sloan Digital 
Sky Survey (SDSS Aba zajian et al.||2009|, the Baryonic 
acous tic (Oscillation Sky Purve y (BUb!S Schlegel et al.. 
2007| and the EUCLID survey ( [Laureijs eral.||201l| . in' 
Section |6j we conclude. 

2. COSMOLOGY WITH VOIDS 

In this Section, we recall the basic equation at the base 
of the Alcock-Pasczyhski test ( Alcock fc Paczynski|1979 1 
applied on voids. This test conies comes from the relation 
between the comoving angular distance Da and the red- 
shift z of an event in a Friedmann-Lemaitre-Robertson- 



Walker (FLRW) cosmology: 
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with rifn, the mean matter density, and Q\ the Dark 
Energy density, both normalized to the present critical 
density. For sufRcientlysmall curvatures fc , it is possible 
to invert the relation Q to derive the redshift from the 
comoving distance r 
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If we look at a cosmic object, e.g. a galaxy, or a cosmic 
void or the BAO, at redshift z, it has an extent Sz in the 
redshift direction and 5r in the angular direction defined 
as 

Sr = DAiz)Se. (6) 

Additionally 6z corresponds to a comoving distance 
along the line-of-sight given by the differentiation of 
Eq. M: 
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with the first derivative of f^. We have indicated in 
the introduction that we assume the Universe is locally 
isotropic. Consequently, the large-scale structures must 
not have a preferred direction in average. If we con- 
sider a void consisting in an infinite average number of 
stacked voids of a specific volume, this "stacked void" 
should have the same extent in all directions. We can 
thus assume that 61 = 6r. This equality yields 

Sr _ 

Sz ~ H^E{z)- 



(8) 



which in terms of the projected separation 5d =^ cz56/Hq 
gives 



Sd 



Ho 

c 



Da{z)E{z) ^ Ho 
zflixiz)) 



r{z)- 
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We propose observing this quantity through measuring 
the shape of stacked voids in redshift space as a function 
of redshift. Note that this observable depends directly 
on E{z) rather than through an integral, as is the case 
for methods based on angular diameter distance, probed 
by observing the angular scale of the BAO peak as a 
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Figure 1. Storing voids in a hierarchical tree. The spatially adjacent sub-voids are assembled in a single void, the "parent void", which 
contain at least one additional basin. The identity of the parent void is inherited by the "sub-void" with the smallest core density. The 
parent void is itself part of a "greater void" . The tree is ordered in scale. 



function of redshift, or the luminosity distance which is 
probed by supernova surveys. We expect this to enhance 
the sensitivity of our method to the physical properties 
of dark energy, which appear directly in E(z). This re- 
lation is not strictly an image of the Hubble constant at 
different z as it is modified by DA{z)/{zf'i,{x{z))), which 
is close to one at low redshift. However, it is a good 
proxy for it and it is possible to obtain the equivalent of 
an Hubb le diagram fo r voids. Eq. ^ was already de- 
rived by Ryden ( 1995 ) for universes with no curvature. 
The measurement of the isotropy would be a clean way 
to measure finely the cosmic expansion. 

3. FINDING AND STACKING VOIDS 

In this Section, we present the algorithm that we de- 
veloped to locate voids and stack voids on an expanding 
metric. In the following, we will use the following con- 
vention. The effective radius of a void corresponds to the 
radius of the sphere of equivalent volume. So, if V is the 
volume of the void, 



redshift Z, we define: 



in 



-V 



1/3 



(10) 



First, in Section |3.1[ we define the coordinate system 



that we use in this work. In Section 3.2, we give the 



details of the algorithm to find and choose the voids that 



are of interest for the stacking procedure. In Section 3.3 
we describe the stacking procedure. 

3.1. The coordinate system 

The fundamental point on which we base our method 
is the capability to find and stack void structures, even 
if they are strongly distorted. For simplicity, we adopt 
the following definition, which a posteriori we will show 
is robust to distortions in the coordinate system. We 
will use the infinite remote observer approximation for 
the redshift coordinate and the planar approximation for 
the angular coordinate. We consider density tracers, e.g. 
galaxies, in a hybrid coordinate system {x,y,z). For a 
tracer t, e.g. a galaxy, at the sky position (9, 0), 9 being 
the sky latitude and (j) the sky longitudes, and located at 



cZ 

Ho 
cZ 

H~o 
cZ 



cos((/)) cos(6') 
sin((/)) cos(6') 



sin(6') 



(11) 
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In the infinite remote observer approximation 9 ^ t:/2 
and thus z ^ j^, but other directions are adequate pro- 
vided that the extent on the sky is small. 

3.2. Organizing voids in tree 

From the volume sampled by the tracers {t}, we ex- 
tract a box of side L in the (x, y, z) coordinate. We 



now use the ZOBOV (Neyrinck 2008) algorithm to com- 
pute and locate local minima m the density of tracers, 
assuming they are a sample of the underlying matter den- 
sity field. ZOBOV finds local density minima on density 
field sampled by particles and their associate catchment 
basins, in the lang uage of the watershed transform (e.g. 
Platen et al. 20071. It does so using a Voronoi tessella- 
tion, derived from the Delaunay tessellation applied on a 
set of tracers. In addition, basins are assembled in voids, 
starting from the lowest density basin, such that: 

- each basin is a void 

- two basins are assembled in one void if they share 
a common boundary, and that the density on this 
boundary is the lowest for each of the void. Basins 
are always assigned to the voids which have the 
lowest core density. 

The whole volume sampled by mass tracers is thus par- 
titioned into a set of basins. Each basin corresponds to a 
void which itself is a collection of basins. Thus, the voids 
naturally acquire a tree-like structure for which voids 
have both a single "parent" and possibly many "chil- 
dren" . We define an order on the tree such that a void 
is the immediate parent of another one if it shares the 
same zones as the child void and at least one more. We 
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Figure 2. Void tree from an N-body simulation. We show the 
convex hulls of particles within a AOh~^ Mpc deep slice, and be- 
longing to a void. A void is visualized if it has either a mean 
relative density less than —0.4, for clarity of the representation, 
and an effective radius within [8; 50]/i~^ Mpc, or is within a sub- 
tree of an drawn void. The color encodes the depth of the void in 
the tree as indicated by the top color bar. 



no te that this corresponds to a different tree than defined 
by Aragon-Calvo et al. (2010). Also, contrary to many 
galaxy based void hnder, tins method does not a ssume 
a hard-coded void shape (see Colberg et al.| |2008[ for a 
review of void finders). 



It does not rely on a smoothed 
density field but only on the topology of the tracers, such 
as particles in A^-body simulations or galaxies in observa- 
tions. This is not the sole technique. More recently the 
complete topology of the cosmic web which can be de- 
riv ed through a s et of tracers has been formally studied 
by Sousbie (2011 1 in the context of the persistence. The 
tree that we have developed helps separating overlapping 
voids such that the same volume is not used multiple 
times in a statistical analysis. 

In our approach, there may not be a natural single root 
for this tree because we do not have periodic boundary 
for the box side L. So we introduce an artificial root 
node for which all child nodes correspond to the parent- 
less voids. Additionally, we compute the mean density of 
each void in the tree, which is a non-monotonous function 
of the depth in the tree. 

3.3. Void stacking 
Voids have complicated shapes mainly produced by the 



gravitational s hear field (Park & Lee 2007 Lavaux & 
Wandelt 20101. It is therefore ditticuit to use them in- 
dividually as a probe of the effect of cosmological ex- 
pansion. However, using the assumption of the isotropy 
of density fiuctuations, we expect that the average void 
shape is spherical in physical coordinates. To achieve 
a correct stacking, at least three important features are 
required: 



1. for a given stack, the algorithm must only select 
voids within a narrow volume range. This condi- 
tion will cause a specific average void shape of the 
given size to emerge from the stacking. 

2. each void of a stack must be unique and must not 
overlap with another void. Enforcing this condition 
allows the removal of spurious correlation in the 
void shapes, which could systematically affect the 
result. 

3. each void in a stack must be centered precisely on 
its average lowest density. This is necessary to min- 
imize the effect of halos, which are expected on the 
boundaries of voids. These halos may bias the po- 
sition of the center reducing the amount of signal 
available from the average void shape. 

The first requirement is satisfied by the ZOBOV void 
finder. ZOBOV provides a volume for the void which is 
exact from the point of view of topology. The second 
point is achieved by stori ng th e voids in the tree struc- 
ture mentioned at Section [3.21 

The center of the void could be derived directly from 
the Voronoi tessellation. However we expect this inferred 
center to be unstable with respect to shot noise in obser- 
vations, or the Lagrangian grid in A-body simulations. 
These would cause spurious effects in the stacking pro- 
cedure. We opt for computing the mean lowest density 
position using the volume- weighed barycenter of the trac- 
ers attached to a single void of the stack. The Voronoi 
tessellation gives a volume Vi surrounding each tracer i. 
We compute the average position xy of the center of the 
void V by computing 



1 



(14) 



where i runs over the tracer, e.g. galaxies, in the void 
V, and Xj it s position in the coordinate system given in 
Section 13.11 

Finally, some voids of the stack have large clusters near 
their estimated center. In the limit of an infinite number 
of voids, we expect these cases to be vanishingly small. 
However, we do not have an infinite number of voids in 
reality. We have decided to avoid these cases by enforc- 
ing that the core density of the void should be really 
empty. We define the core density as the mean matter 
density within a sphere of fiducial radius rcff/4. This 
comes at the cost of a smaller numb er o f available voids 
in the stack. A posteriori, in Section [4Tj we see that our 
adopted fiducial radius is smaller than the actual size of 
the "core" of the stacked void. 

To summarize, we proceed as follows: 

1. We put t he tracers in the coordinate system of Sec- 
tion [O 



2. We extract a parallelepiped volume of side Lx Lx 
Lz which resides entirely within the region spanned 
by the tracers. 

3. We normalize the coordinates to 1 in each direction. 

4. We apply the Zobov void finder algorithm. 
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Figure 3. Void stacking in the final, z = state of the simulation. We show the result of the stacking procedure given in Section |3.3| 
when it is applied to all voids of size 8 Mpc in the full volume of the simulation without introducing cosmological and peculiar velocities 
distortions. Left panel: a filled co ntour plot of the density n{d,z). Contours above the density 1.2 are shown in light gray. We fitted the 
simple ellipsoidal model of Section |4.3| on the binned density to estimate the ratio between the axis along the XY direction (x-axis) and 
the axis along the Z direction (y-axisj. Right panel: the three-dimensional average density profile, in shells, of the stacked void (solid line) 
and the fit to a cubic density profile (dashed line). The horizontal solid red line is the mean density. 




Figure 4. Stacked void profiles at z = - We show the profiles for voids selected with different effective radii. Left panel: mass density 
profiles in thin shells. Right panel: average mass density profile in the sphere of the given radius. 
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5. We sto re th e found voids in a tree according to 
Section [X2] 



6. We walk the tree, starting from the root and stop 
whenever the effective radii, r^s, is within a given 
range, between iJmin and i?max- i"oS is obtained by 
computing the radius of the sphere which has the 
same volume as the void. 

7. We compute the position of the particles according 
to the volume- weighed barycenter of Eq. (14) for 
each void. 

8. We compute the density of the void within a sphere 
of radius rcff/4. We only accept voids that have 
a core density less than 20% of the mean matter 
density of the universe. 

9. Around each volume-weighed barycenter, we ex- 
tract a spherical volume, in the coordinates of Sec- 
tion [XT] of radius i?cut = 3 x Rmax- The center of 
the extracted volume is put at the origin. We do 
this for all selected voids. 

The resulting particle distribution gives what we call a 
stacked void. We transform coordinates from (x,y,z) to 

{d — y^a^+'j/^, \z\). We bin this distribution and divide 
the resulting density by the number of stacked voids, 
the bin width and d, the Jacobian of the transformation. 
This procedure yields the density n{d,z) of tracers per 
unit volume per void. An example of such a distribution 
is given in Figure [3] 

4. TESTS ON iV-BODY SIMULATION 

We have run a series of three pure dark matter iV-body 
simulations with different realizations of the initial condi- 
tions but for the same cosmology. The dark matter parti- 
cles have different sampling properties than the galaxies. 
In particular, galaxies are biased tracers of the matter 
density field. As we are relying on topological proper- 
ties of the large-scale structures, the result should be the 
same for the two. We discuss in Section|6]thc limits of our 
approach. The volume of each simulation is given by a 
cube of side L = bOOh^^ Mpc. Each simulation has N — 
512^ particles. We have adopted a ACDM-WMAP7 cos- 
mology with the following parameters: fJb/i^ = 0.02258, 
n^h^ = 0.1108, H = 71 km s^^ Mpc-\ w = -1, rig 1, 
As = 2.34 X 10-9. This corresponds to = 0.045, 
f^M = 0.264, CTg = 0.84. Each particle has a mass 
mp = 2.0510^^/1"^ M©. The transfer function for den- 
sity fluct uations for t his cosmology is computed using 
CAME (iLewis et al.|l2000). The initial conditions are 



generated using ICgen[] which uses the transfer func- 
tion to generate a density field from the primordial power 
spectrum. 

4.1. Stacking voids in pure comoving coordinates 

In this Section, we consider the ideal case of voids 
stacked in comoving coordinates,! e. we purely consider 
the distribution of dark matter particles as given by the 
A^-body simulation. We consider voids for which roff is 
between iJmin — 8h~^ Mpc and i?max — 9h^^ Mpc. 



We give the result of the void stacking procedure 
for one of the A^-body sample realizations in Figure [3j 
In the left panel, we show the density profile of the 
stacked void, where on the horizontal axis corresponds 
to d = \/ + y^, with X and y the first and second co- 
ordinate of the particle in the stacking, and the vertical 
axis to |z|, the third coordinate of the particle. The solid 
black contour corresponds to the res ult g iven by the like- 
lihood analysis described in Section 4.3 The filled color 
contours have been chosen equi-spaced from pjp = to 
pIp = 1.1. The solid line in the right panel shows the 
three-dimensional density profile, in thin shells, of the 
stacked void. 

The one-dimensional density profile shown in the right 
pan el of Figure [3| is s i milar to the one shown in Figure 



2 of Ceccarelli et al. (20061, though we use pure dark- 
matter simulation in place of mock/real galaxy samples. 
The matter shell around the void is clearly visible for 
radii greater than the chosen i?niin- 

The inspection of both panels shows a number of in- 
teresting features. It is clear that the voids stack coher- 
ently and form a region of lower density for \J <P- ^ z'^ < 
8h~^ Mpc. The mass density inside the stacked void, in 
both panels, is featureless. Outside the expected bound- 
ary of the stacked void, at V d'^ + z"^ > 8h~^ Mpc, the 
density profile continues upward to p ~ l.lp and then 
falls back to homogeneity. This shell is clearly seen in the 
one-dimensional averaging of the density profile shown in 
the right panel of the same Figure. Our d ensity profile 
has sim ilar features to the on es p resented by|Be nson et al.| 
Padilla et al.l (|2005|) and I Ceccarelli eTal. (2006^ 



n the two-dimensional mass density diagram ot the left 
panel of Figure [3} it is clear that this is a real spheri- 
cal shell, and that it harbors clumps highlighted by the 
white regions at p/p > 1.3. 

Our profiles are not as st eep as in observations and 



some mock galaxy catalo gs (Benso n et al. 2003 Padilla 
et al.||2005[ |Hoyle fc Vogeley, 2004| e.g. J. I'his may be 
due to our use ot pure dark-matter s imulation in place of 
mock galaxy samples. The results of Benson et al. ( 2003 1 
(Figure 11) also indicates that observing tamter galaxies 
tends to smooth the density profiles, as expected. Even 
though for the purpose of applying to observations (such 
as in |Ceccareili et al.||2006| it may be better to use their 
form, we have empirically found that a cubic function fits 
the simulated density profiles adequately for the purposes 
of this analysis. We use the following form 



p{r) 



= Ao + A3 



Rv) 



(15) 



Available from littp://www. iap.fr/users/lavaux/ 



with i?v the radius of the stacked void, Ai the parameters 
against which the fit is computed. As we have voids 
with 8h~^ Mpc< Tcff < 9h~^ Mpc, we have set Ry = 
8h-^ Mpc. The best fit gives Aa = 0.13 ± 0.01, A3 = 
0.70 ± 0.03 for voids of 8h-^ Mpc. 

First, we note that the value of Aq attributes a signif- 
icantly non-zero density to the center of the void. This 
may be due to a resolution effect as putting one par- 
ticle within (l/j-^Mpc)'^ yields a density fluctuation of 
S = —0.2 for the resolution of our simulation. On the 
other hand, this value is the same for large voids where 
resolution effects should be milder. This was already 
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noted by [Colberg et aL] (|2005 1 . 

Second, we note that (1 — Ao)/A3 is not strictly equal 
to one. This can probably be explained by the size of the 
interval of accepted effective radii. We exp ect that large 
bin sizes, such as in Colberg et al. (2005), may signifi- 
cantly change the prohie ot the stacked void, notably at 
r ~ Ry owing to the non-commutativity of the operation 
of re-scaling and averaging. This problem also affects all 
the quantities that we may derive from this profile, like 
peculiar velocities. The result of the fit (dashed line) and 
the actual void profile (solid line) is shown in the right 
panel of Figure [31 Visual inspection shows a good agree- 
ment between those two profiles within the bounds of the 
voids. We will use this profile for the remainder of this 
work. For our method the choice of the profile is only 
important insofar as it does not bias (or remove signal 
to noise from) the shape measurement. 

In Figure |4] we show the density profiles, both in shells 
and cumulative, with radii normalized to the sought ef- 
fective radii of the voids in each stack. We have consid- 
ered voids with roff between 4— 5/i~^ Mpc, 8 — 9/i~^ Mpc, 
10- ll/i-i Mpc and 14- IS/i"^ Mpc. We clearly see that 
all voids with rcff greater than 8h~^ Mpc have nearly ex- 
actly the same density profile, even at radii larger than 
i?v- While th e universality of void profiles was previ- 
ously noted by Colberg et al. ( 2005 ) , they did not show 
that this universality also extends outside the void, in 
the shell region. Our density profile s compares well also 
with that in the seminal paper by [van de Weygaert fc| 
van Kampen p993, ). Note that their void size definition 
is radically difierent from ours. Finally, from the left 
panel of Figure |4j we note the presence of a structure 
looking like a "core", which does not extend above half 
the effective radii of the individual voids. This justify a 
posteriori our choice for the mean density within roff/4 
to be sufficiently low to accept a void in the stack. 

4.2. Stacking with redshift distortion 

In this section we consider the possibility of applying 
the algorithm for void detection of Section [3] We have 
selected five snapshots for each of the three simulations. 
They correspond to the simulated universes at = 0, 
Zb = 0.25, Zh = 0.54, Zb = 0.81 and Zb = 1.0. We have 
considered that the first two dimensions of the positions 
of the particles in the snapshots correspond to the an- 
gular coordinates and the last, zq, to the distance along 
the line of sight, relative to the face with zc = 0. As 
our simulated universes are flat, we have converted the 
comoving positions zc in redshift positions z using 



z = x 



— {zc + x{z^)))+-, (16) 
c J c 



with X as defined in Eq. ([3|, the peculiar velocity of 
the particles in the z direction, c the speed of fight. We 
have not included the additional Da{z)/z term present 
in e„(z) (Eq. [9| which adds an additional small effect 
at low redshiftT^ In each case, we have extracted a box 
of AQOh-^ Mpcx400/i~i Mpcx40,000 km s^^ from the 
distribution of particles in redshift coordinates. We have 
run the void identification and stacking algorithm on the 
particles of this box. Additionally, we have considered 
the stacked particles when peculiar velocities are either 
included (dubbed "mock catalog with full redshifts") or 



excluded (dubbed "mock catalog with pure expansion 
redshifts"). 

In Figure [5j we show the result of the stacking for one 
of the simulations, for the snapshot at redshift z = 1. 
For visualization purposes, the effects of expansions were 
removed after having stacked the voids. In the left panel, 
we show the result of the stacking algorithm for the mock 
catalog with pure expansion redshifts. In the right panel, 
we show the same test but for mock catalogs with full 
redshifts. We highlight with black contours the fitted 
three-dimensional density profile of Section [4.1[ The plot 
is corrected for the inferred expansion, so the contours 
are perfectly circular. Qualitatively, they match the color 
coded density in the left panel. 

We note in the right panel of Figure [5] that there is 
a non-trivial deformation of the void. Atlow (d, z), the 
void is emptier and slightly elongated in the redshift di- 
rection. This is expected because of a finger-of-god effect 
in the void. At high (rf, z), the void is flattened. This 
is clearly seen by considering, e.g., the yellow iso-density 
in the right panel. The decrease of the peculiar veloc- 
ities are not sufficient to explain the amplitude of this 
effect. We have tested our fitting procedure on mock 
voids whose shaped have been transformed by the aver- 
age peculiar velocity field both expected and measured 
in voids. We have found that they are not introducing a 
significant pancaking effect. 

We have investigated the origin of this systematic ef- 
fect, which happens to be non-trivial. It may be ex- 
plained by a two-step process. First, large-scale redshift 
space distortions induce a selection bias on cosmic voids: 
cosmic voids with collapsed structures along the angular 
coordinate direction are slightly disfavoured when they 
are binned by void sizes. This discrimination is generated 
by large scale flows that induces a small modification of 
the void size, which is sufficient to displace voids from 
one void size bin to another. However, because voids do 
not have a flat size distribution, it is statistically much 
less likely to displace a large void into the bin than to 
displace it out to a bin corresponding to larger voids. 
The stretching is far less present when large clusters are 
present along the line of sight, which means that these 
voids stays in the same void size bin. This selection ef- 
fect makes the distribution of large-scale structure out- 
side the void slightly anisotropic. Additionally, as halos 
are preferentially located at small angular distance and 
high redshift distance from the void center, the finger- 
of-gods that they produce cause a thickening of the void 
wall in the redshift direction, which in turns cause the 
pancaking of the voids. 

At smaller radii, like ~5 h^^ Mpc, the stacked void 
is seemingly spherical because it is not contaminated by 
any of the two above effects. By considering the void 
shape near the half-radius from the center therefore min- 
imizes possible biases due to peculiar velocities. This is a 
robust procedure, but it is not lossless. Further modeling 
of the profile could improve signal to noise in our shape 
inference, described in the next section. In Section 5.1 
we propose a simple alternative de-biasing scheme. 

4.3. Void shape inference in redshift space 
Stacked voids have an isotropic cubic density profile in 



comoving coordinates, as found in Section 4.1 We model 
the redshift space distortion by fitting the density n(d, z) 
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Figure 5. Impact of peculiar velocities. Density in stacked voids in redsiiift coordinates, in the snapshot at 2 = 1, for a 
400/i^^ Mpcx400/i~^ Mpc slice with a thickness of 10,000 km s~^. Left panel: expansion only, without the peculiar velocities in the 
redshift positions of particles. Right panel: we have introduced the distortion due to peculiar velocities. The black contours are the same 
in the two panels and highlight the fitted density profile to the left, expansion-only case. 



estimated using the stacking procedure of Section [3. 3[ in 
redshift/angular coordinates this time, to the function 



n(d, z) 



fiQ^ {{d/adf + {zla,ff''^ ,71- 



with riQ the density at the minimum in (c?, z) = (0,0), 
ad the semi-axis along the angular coordinate direction, 
ttz the semi-axis the redshift direction, rimax a maximum 
density value. We expect rimax ought to be near unity 
to show convergence to the mean density. However, we 
leave it as a free parameter because of the presence of the 
shell around the void and the limited accessible volume 
around the stacked void which could bias this value. 

To fit the model to the estimated density we assume 
that the fluctuations according to the ellipsoidal model 
are Gaussian but with two different variances depending 
on the location in the void. The likelihood takes thus 
the following shape: 



X^("-o, "-max, ad, a^,aa, cti) 

i=i j=i 



{n{di,Zj) - n^^jY 
a'^{di,Zi) 



+ 2 \og{(7 {d,, Zi)) j , 
(18) 



with the z-th and j-th bin, for which the {d, z) take 
the value (d^, Zj), riij the value estimated in the bin 
Nd the number of bins in the d direction, the 
number of bin in the z direction. Sij is either one or 
zero depending whether we want to include the bin (i,j) 
in the optimization. We limit the fit to the disc of radius 
Rcut, as there is no stacked void data at distance larger 



than this. Sij takes the form 



1 if ^d1 + {z,/EY <R, 
otherwise 



(19) 



where we may correct for the expansion in the pixel se- 
lection through the coefficient E. In practice, we keep 
E' = 1 in all the following. 

At the position {d,z), and within the void, we enforce 
that the Gaussian part, a, of the distribution 7?. scales as 
1/^/d. This follows from the cylindrical averaging when 
building the stack of voids. Outside the void, we take 
a fixed standard deviation to account for uncertainty in 
approximating the outside profile by a single constant. 
The standard deviation is spatially varying according to: 



a{d, z) = 




ih- 



^ iin{d,z) < 



otherwise 



(20) 



We find the parameters (no, "-maxi Od, ctqi i^i) and their 
error bars by running a Monte-Carlo Markov-Chain ex- 
ploration of the four parameters on the sub-sample of 
pixels for which ^J~^^Y^r\zjJ~EY < Rent ■ AH the measure- 
ments of ratios shown in the Section [5] are made using 
this technique. 

The error model contains two components — the Pois- 
son error due to the number of tracers in each pixel, and a 
correlated error due to dense clumps which occurs around 
individual voids and which have not been completely re- 
moved by the averaging procedure. Such clumps of scale 
< lh~^ Mpc are visible in Figure [s] A crude way to 
model such correlated fluctuations is to choose the pixel 
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Figure 6. Result of the likelihood analysis on the binned parti- 
cles of the stacked void - We give in solid b lack line the iso-density 
contour of the stretched cubic model (Eq. |17| fit using the likeli- 
hood analysis of Section |4.3| For this Figure, we have considered 
voids with r^ff = 8h~^ Mpc. We plot the underlying density field 
between the null density and l.lp, with p the mean density of the 
slice. Density pixels are 2h~^ Mpc. The error is modeled as inde- 
pendent from pixel to pixel. 

size to be large enough that such fluctuations mostly af- 
fect single pixels only. We conservatively choose a pixel 
size of 2h~^ Mpc to satisfy this criterion. For pixels of 
this size the Poisson error due to individual tracers is 
entirely negligible. 

We illustrate in Figure [6] the outcome of the likelihood 
analysis on a stretched stacked void, obtained from a 
mock light-cone at z = 1. We show both the density 
field of the stacked void and the iso-density contours of 
the fitted profile using the likelihood analysis. The iso- 
density contours follow the outer edge of the void. The 
fluctuations of the binned density field looks clearly ran- 
dom and uncorrelated at this resolution. This was not 
the case in Figure [5j We conclude that the likelihood 
thus behaves as designed. 

We tested the robustness of this procedure to changes 
in the number density of tracers by re-running the en- 
tire pipeline on a subsample of the A^-body particles. 
The results essentially did not change if the number of 
particles was reduced by a factor 5 (to ~ 0.2/i'^Mpc~'^), 
which corresponds to a typical galaxy density expected 
for the EUCLID survey at low redshift (z < 0.1). This 
confirms that Poisson error due to the number of tracers 
is negligible for the large pixel sizes required to reduce 
pixel-to-pixel correlations due to non-linear structures. 

4.4. Number of voids 

We give in Figure [7] the averaged number density of 
voids, in comoving coordinates, for the three simulations, 
at each redshift and each Toff that we have considered for 
estimating the Hubble constant. 
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Figure 7. Comoving number density of voids. We show results 
for purely expanding universes, the result when peculiar velocities 
contaminates redshifts is the same. The dashed lines give the lowest 
and the highest value of the density when the three simulations are 
considered. 

The dependence of the number of voids with redshift 
is as expected. Small voids should be more abundant 
at high re dshifts and large voids more abun dant at low 

The rela- 
As noted 



redshifts (Sheth & van de Weygaert 



20041. 



tion pivots about a radius of ~8/i~^ Mpc. 
by these same authors, the number density of voids de- 
pends principally on their volume and then on redshifts 
through growth of structures. The number of voids de- 
tected, when peculiar velocities contamination is added 
in our mock catalogs, is roughly the same and with the 
same dependence with radii. We note a small but sys- 
tematic destruction of voids of 4:h~^ Mpc. It is plausible 
that the original topology is lost at these small scales be- 
cause of the contamination by fingers of god, which can 
be effectively as deep as 10h~^ Mpc. 

In FigurefT) we note that whatever the dependence of 
the physicalnumber density of voids with redshifts, for 
voids as we define them in redshift space, this dependence 
is small within the redshift range z = 0—1. Thus for the 
purpose of this work, we neglect the time dependence of 
the void abundances and focus on the scale dependence. 
For each radius, we average the densities at all redshifts 
and fit an exponential law, which seemed most suited for 
representing this set of curves, on the number density as 
a function of res, for the range 4 — 14h^^ Mpc: 



1 /i3 Mpc"^ ~ 
(3.5 ±0.2)10^3 



-(0.632 ± 0.006)—^^ 



1/i-iMpc 



(21) 

We show this relation in Figure |7] with a thick black 
solid line. This relation is used for the Fisher-Matrix 



analysis of the Section 5.2 as an approximation of the 
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Figure 8. Hubble diagram derived from the voids. Hubble diagrams derived from voids witli three effective radii, without including 
distortions due to peculiar velocities: 6/t~^ Mpc (top- left), 8h~^ Mpc (top-right), 14/i~^ Mpc (bottom- left). The bottom-right panel 
corresponds to rgft = 8h~^ Mpc and including distortions due to peculiar velocities. We show the actual expansion in the mock catalogs 
(black line) and the recovered average expansion from stacked void shapes (colored dashed lines) for the three Af-body simulations. The 
colored error bars show the standard deviation inferred using our statistical model described in Eq. |18| In the bottom right panel, we used 
a 40,000km s'^ thick slice instead of 10,000km s~^. 
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behavior of the number of voids as a function of scale. 
We estabUsh this empirical relation without fitting to any 
void formation mod els, such as in the one in'Sheth & van 
de Weygaert ( 2004 1 . This relation only reflects how voids 
are detined by our algorithm described in Section [3] 

5. ESTIMATING THE EXPANSION HISTORY 
USING VOIDS 

With what precision can we obtain e^(z), which is 
closely related to H{z), by only considering the redshift 
shape deformation of stacked voids? We use the algo- 
rithms and results of Section |3] and |4] In Section 15. 1[ 
we present and discus s the quality of the derived Hubble 
diagrams. In Section |5.2[ we construct a Fisher-matrix 
analysis of the constraints that can be obtained on Dark 
Energy with this method, and compare it to expected 
baryonic acoustic os cillations constraints from the Dark- 
Energy Task Force ( [Albrecht et~al]|2006 ). 



5.1. Results from simulations 
We show in Figure |8] the complete "Hubble" diagram 



obtained from voids ofradius 6h~^ Mpc (top-left panel), 
8h~^ Mpc (top-right panel) and 14/i"^ Mpc (bottom- left 
panel). These voids were obtained from the simulations 
presented in Section [4] using the algorithm described in 
Section [Sj For all these diagrams, we have not contam- 
inated tne cosmological redshifts by peculiar velocities. 
Thus they simulate the measurements in a purely ex- 
panding universe. We show in the bottom right panel 
the Hubble diagram for voids of 8h~^ Mpc, when pecu- 
liar velocities are included in redshifts. We show the ac- 
tual measurements on the three A^-body samples, along- 
side with the 68% error bar as derived from the Bayesian 
redshift shape inference. Each color corresponds to the 
same simulation in all panels. 

We have used five snapshots of the simulation, cor- 
responding to time, expressed in redshift, z = 0, z = 
0.25, z = 0.538, z = 0.818 and z = 1. For each of 
these snapshots, we have appli ed cosmological expan- 
sion as explained in Section 4.2 For voids of 6h~^ Mpc 
and 8/1-1 y^^^^ 

we divide the volume in four slices 
of 10,000 km according to the redshift direction. 
For voids of 14h~^ Mpc, we keep the full volume of 
400h'^ Mpcx400/i-i Mpcx40,000 km s~^. 
_ In black solid line, we plot the average local stretching 
E{z,Az), derived from Eq. (|8|, 

fZ+Az 

(22) 



1 

W) 



1 



dz 

wr 



where z is the minimum redshift of the slice and Az its 
thickness. Assuming that voids are uniformly distributed 
in the volume, the stacked voids should be stretched by 
E. The theoretical expectation and the measurement in 
A^-body simulation are in agreement for the three void 
sizes Qh"^ Mpc, Sh'^ Mpc and lAh'^ Mpc, as clearly 
shown by the comparison in Figure [s] (black solid line) . 
We note that the standard deviation derived from the 
posterior are quite larger than the residual for the voids 
with Toff = &h~^ Mpc. This is linked to our binning 
choice and the need for a better model of variations of 
the density at scales smaller than 2h~^ Mpc. 

Considering the Hubble diagram obtained when the 
distortion due to peculiar velocities are included, we note 



1.8 



Simulation R=8 h ' Mpc 
Tlieory 



0.6 
Redshift 2 



Figure 9. De-biased Hubble diagram - Same as Figure [8] 
We siiow tile Hubbie diagram derived from voids v^ith r^ff = 
8h~^ Mpc, including peculiar velocities distortions and after cor- 
rection of the systematic bias. 



an important systematic effect in the bottom right panel 
of Figure [81 This distortion is a consequence of the pan- 
caking effect that has been discussed in Section 4.2 It is 
unfortunately non-trivial to model. 

The pancaking effect does not affect the overall stretch- 
ing but it transforms the structure of the density field 
outside the void, which may bias our statistical estimator 
of the shape. It is also expected to be weakly dependent 
on cosmology, essentially cts and fim which both affects 
the amplitude of peculiar velocities inside and outside 
clusters. We have estimated using mock catalogs from 
which finger-of-gods were removed that it corresponds 
to a ~5% effect on the stretching. The thickening of the 
walls by the finger-of-god corresponds to the rest of the 
apparent flattening (~11%). On one hand, this system- 
atic should depend slightly on voids sizes as larger voids 
have a comparatively smaller flattening induced by the 
finger-of-god. On the other hand, because of the much 
smaller number of voids, the error bars looks sufficiently 
big to hide this effect in the present study. We propose 
to address this issue in a future work focused on mock 
galaxy catalogs. 

As a first approximation we take the bias induced by 
the pancaking as constant, and leave more detailed sta- 
tistical and dynamical modeling for future work. We 
empirically find that the observed stretching E{z) (ci, in 
observations) should be multiplied by the constant de- 
biasing factor of 1.16 ± 0.04. We have checked that this 
constant is not strongly dependent on void radii. We 
note that this multiplicative factor is unimportant for 
the actual determination of the equation of state of Dark 
Energy. This statement is correct as long as the pancak- 
ing effect is not redshift dependent, which could alter the 
slope of the relative between redshift and the stretching. 
Consequently, though we adopt a factor of 1.16 ±0.04 for 
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this section, this factor should be left as a free parameter 
in any attempt to fit the observations of void eUipticities. 

We show the resulting diagram in Figure |9] We note 
that the effect has now disappeared except at the lowest 
redshift where this leads to a slight overstretching. This 
overstretching is due to a competing effect only present 
at small redshift and for small voids. 

5.2. Fisher-matrix analysis for Dark Energy 

To assess the power of the Alcock-Paczyhski test using 
stacked voids, we derive the Fisher matrix for the Dark 
energy properties that can be derived from the Hub- 
ble constant. We consi der a Chevalier-Polarski-L i nder 
(GPL) pa rametrization ( Ghevallier fc Polarski|200T Lin- 
der||2003 ) of Dark energy equation ot state: 



w{z) — wo + ' 



'l + z 



(23) 



For sufficiently low redshift the reduced Hubble constant 
is thus 



E{z,wa,Wa) 



(1 + zf + Ha exp 



-3 / dz 
'o 



, l + w{z') 
l + z' 



1/2 



(24) 
We 



with w{z) as defined in Eq. (23), and oj,,! = ^mh^ ■ 
do not assume that the cosmology is flat. 

We consider a hypothetical measurement of the 
stretching constant ev(z) from the study of voids. Thus 
we have an estimate of ev(z) at different redshifts zi that 
we label ev.i- Each of the estimates has an independent 
random error variance Vi. The likelihood C of the cos- 
mological parameters p is thus simply described by: 



(25) 



In p, we include a parameter, which corresponds to 
the overall flattening of the voids due to effects of redshift 
space distortions. We use the value and the uncertainty 
as determined in Section [5j This prior is of no conse- 
quence on the Fisher-Matrix derived for the EUGLID 
survey as it is sufficiently dense to constrain its value at 
lower redshift. The definition of the Fisher matrix is 



dC dC 
dpk dpi 



(26) 



with k,l € p and the averaging is taken over all the pos- 
sible realizations of the data-sets {ev,i} given the cosmo- 
logical parameters p and the noise determined by {Vi}. 
In our case, this definition simplifies to: 



F, 



k,l 



E 



1 dev dcv 
Vi dpk dpi ' 



(27) 



To complete the evaluation it is required to have an 
estimate oi Vi = {{e^s — ev(^i))^), the variance of the 
estimated Hubble constant in each redshift slice. This 
variance may have the following typical dependence: 



V,^ \Y,mj,N{R„z))f 



(28) 



with e{R, N) the standard deviation of the estimator if 
we stack TV voids with a size R, n{R, z) the number den- 
sity of voids of size R at redshift z. In Eq. (28), we are 



summing over all the void effective radii that are observ- 
able in the slice i. If we bin void by size, they may form 
statistically independent stacks. It is thus possible to im- 
prove the variance of the local expansion factor E{z) by 
considering all possible sizes at once. From the tests that 
we have run in Section |4j we know that binning voids in 
bins of A = lh~^ Mpc in effective radii gives adequate 
results. Thus, in Eq. (pSl), we use 



i?7 — i?n 



ohs.{z) + jA, 



(29) 



with j < (i?„iax obs. - Rmin obs.)/A. Of coursc, it is not 
possible to observe void which have a size larger than 
half the width of the slice i. Gonsequently, 



R 



max obs. 



c6z 
2Ho' 



(30) 



with dz the width of the slice, expressed in redshift units. 
We note that this is only an approximation. 

We have chosen the following scaling fo r th e minimal 
observable void size i?inin obs. used in Eq. (129) 



R, 



mm obs 



(z) = min (eh^^ Mpc; Srn(z,)"^/^) , (31) 



with n(z) the mean comoving density of galaxies at red- 
shift z. n(z) typically depends on <f>[M), the luminosity 
function of galaxies in the observed band. For a magni- 
tude limited catalog, as the SDSS Main galaxy sample, 
assuming magnitudes are corrected for evolution n(z) is 
related to (t){M) as 



n(z) 



m-51ogi(,(ci(z)) 



(j){M) dM 



(32) 



Our choice of i?min obs. ensures that the size of the struc- 
tures that are observed is limited by the number density 
of tracers. The Fisher matrix, and of course the fig- 
ure of merit for the determination of Dark E nerg y, is 
going to depend on Sr- We discuss in Section |5.3[ the 
impact of varying Sr- We do not accept voids smaller 
than 6/1""'^ Mpc because they may be easily disrupted 
and strongly contaminated, in redshift coordinates, by 
distortions due to peculiar velocities and are more diffi- 
cult to identify. We leave the exact determination of the 
lower limit of observable voids for future work. 

We approximate the number of voids in the slice of 
thickness 5z at redshift z by 



ky ^Jf^ (^(^ + '^^)^ ^ X{zf) nvoid{R, z). 



47r 
y/s 

(33) 

We determine e, n from the three simulations that we 
have run for a AGDM-WMAP7 cosmology. As this cos- 
mological model gives a g ood fit to GMB and large-scale 
structures observations ( |Komatsu et al.|20TI Reid et al. 
2010), it should give a fair representation of the statisti- 
cal structure of the galaxy redshift survey. 

For the three void sizes that we have considered in this 
work, 6h~^ Mpc, 8h~^ Mpc, and lAh~^ Mpc, we show 



in Figure 10 the function e(i?, TV). We may approximate 
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of sky 
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function 



-20.83 
-1.20 
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M, = ■ 
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I SPSS Collaboration fc Blanton|2000 1 



: 2.63 lO-^" h'^Mpc" 
M, = -19.42 
a = 3.90 



1.7 10^ 



SDSS-DR7 (LRG) 24% 



r = 19.5 0.45 



I Cool et al.|2008 



10^ 



BOSS 



24% 
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r = 20 



0.7 



1.5 10*= 



EUCLID 



36% 



1.1610 
M« = ■ 



-23.39 
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H = 24. 



1.5 



(Kochanek et al.|2001 Jones et al.|2006l 



1.610* 



Table 1 

Survey parameters 



this function by: 



-{R.N) 



0.98 



-0.57 



(34) 



which is independent of the radius R. The reduction in 
the variance is scahng roughly as a Poisson distribution. 
This corresponds to the expectation that the noise in the 
shape estimation comes essentially from the void-to- void 
fluctuations. For future surveys we expect to observe a 
huge number of small voids which would make our esti- 
mate sensitive to the low error cases, corresponding to 
the red points Figure 10 Conservatively, we thus take 



that the standard deviation is scaling exactly as 



e{N) 



1 



N 



(35) 



This slightly overestimates the errors for the shape mea- 
surement and should give a conservative estimate of the 
expected constraints. 

5.3. Application to present and future surveys 

In thi s se ction, we move on to apply the formalism of 
Sectio n 5.2 to three surveys: the SDSS-DR7 (Abazajian 
eF^^OT, the BOS S survey ([S chlcgcl et al. 2007faHd 
the KUClID survey ( Laureijs et al.,,2011) . We use slices 
of 5z = 0.03, which corresponds approximately to a co- 
moving thickness of 'dOh~^ Mpc. This limits the effective 
radius of the void to ^ 45/i~^ Mpc. 

Table [l] lists the parameters of those surveys which are 
useful for the Fisher-Matrix analysis. These parameters 
can be translated, using Eq. ( [32| ), in a number density of 
tracers at a given redshift. This det erm ines the min ima l 
size of observable voids through Eq. ([3l]) . Using Eq. ( 35 ) 



we obtain the number density of voids and with the void 
shape standard-deviation reduction e{R, N). 

For deriving Fisher-matrix constraints, we assume the 
priors that should come out from t he Planck mission . 
We use the Dark Energy Task Force ( Albrecht et al. 2006 1 
prescription for deriving these priors! in addition, we 
apply constraints from Stage II experiments as expressed 
in the DETF report. The figure of merit (FoM) is defined 
as 



FoM 



1 



a{Wa)(j{Wp) 



(36) 



Residual 

Statistical model 
Fit 




Number of voids 

Figure 10. Standard deviation of void apparent shape vs the 
number of s tac ked voids. Shape errors inferred by our statistical 
model Eq. ( |18| (solid lines), the residual between the expected 
shape and the actual best estimate (thick cross markers). We also 
sho w wi th a dashed line the best fit of a power-law on the residual 
(Eq. |34| . The results for the three A'^-body simulations are shown 
here.^ The colors correspond to voids with an effective radii of 
6/i~l Mpc (red), 8/i~^ Mpc (blue) and 12/1" ^ Mpc (green). 



with 
and 



Wp = Wq + {I - ap)Wa 
(SWaSwo) 



1 - Op = 



(37) 
(38) 



Wp is the best estimate we can obtain on the equation of 
state of Dark Energy, evaluated at the scale factor ap of 
our Universe. 

We consider either an experiment consisting of only 
BAO analysis, only Voids analysis or the two combined. 
We do this analysis for the SDSS main galaxy sample, 
the BOSS survey and the EUCLID survey. We give the 
results in Table [2] Figure 11 and Figure [12] 
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BOSS 

SDSS Voids (Main+BOSS), ,i, = 

Both, a, =2 

SDSS Voids (Main+BOSS), ,i, = 
Both », =5 




-0.4 -0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 0.4 

Figure 11. Fisher matrix forecasts for SDSS/BOSS survey. Pre- 
dicted 95% confidence regions, assuming Planck prior and the de- 
termination of tfie Hubble constant H = 72 ± 8 km Mpc~^. 
Tlie Planck priors were obtained using tfie procedure in tfie re- 
port of the DETF. Please note that we plot wg vs Wa and not Wp 
vs Wa- The black solid line gives the constraints derived from the 
analysis of baryonic acoustic oscillations in the BOSS survey. The 
dotted lines gives the constraints derived from the shapes of stacked 
voids using SDSS main galaxy sample and the BOSS sample. The 
dashed line shows the obtained constraint by combining both Bary- 
onic Acoustic Oscillations and void shape measurement. The con- 
straints are either derived assuming Sr = 2 (black) or = 5 (red). 



We note that voids and BAO have roughly the same 
constraimng power when considering the SDSS/BOSS 
experiment (FoM of 71 for BAO only vs. 68 for voids 
only). The combination of the two slightly improves the 
constraints (FoM of 75). On the other hand, voids are 
far superior at constraining the equation of state of Dark 
Energy with the EUCLID survey. We have found a FoM 
of ~1 825 far superior than the constraints derived from 
BAO for the same survey. This comes from the possi- 
bility of using small scales geometry and not only the 
^lOOh^^ Mpc scale which corresponds to BAO. 

The results depend on several parameters that we have 
adopted for the analysis. For example, we have chosen 
.Rmax = 45/i~^ Mpc. We have barely found any change 
in the results by limiting to rcs — i4:h~^ Mpc or by 
changing the thickness 6z of the slice up to 0.1. How- 
ever, changing i?min have a lot more impact on the FoM. 
This indicates that, as expected, most of the information 
comes from the smallest observable voids because of their 
huge abundance relatively to bigger ones. Similarly, we 
have tried varying Sr between 1 and 10 to check the sta- 
bility of the constraints for EUCLID. We have found that 
the Hubble diagram for voids always give a significant 
additional information. As an example, in Figure [T2| we 
show in black (in red respectively) the constraints for 
Sr- = 2 {sr = 5 respectively). The figure of merit may ei- 
ther strongly improve whenever Sr is reduced, up to one 
hundred times improvement, or diminish by a factor of 
a few when Sr is increased. We also tried to check the 
influence of the sky coverage on the FoM. Interestingly, 
reducing the EUCLID survey to 100 deg^ yields a figure 
of merit of ~380 for voids, while the one we expect from 



BAO EUCLID 
Voids EUCLID, ,iv =2 
Voids EUCLID, ,5, =5 




Figure 12. Fis her matrix forecasts for the EUCLID wide survey. 
Same as Fig. [IT] but for the EUCLID survey. The solid (dotted) 
line gives the constraints derived from the analysis of Baryonic 
Acoustic Oscillations (shapes of stacked voids, assuming Sr = 2 in 
black and = 5 in red). 



Method 



Data 



FoM 



BAO 




BOSS 


71 


Voids {sr = 2) 




SDSS-t-BOSS LRG 


69 


Voids {sr = 5) 






68 


BAO-f Voids {sr = 2 


or 5) 


SDSS-I-BOSS 


75 


Voids (sr = 2) 




EUCLID 


~1 825 


Voids {sr = 5) 






~273 


BAO 






98 


BAO-I- Voids {sr = 


2) 




~2 956 


BAO-I- Voids {sr = 


5) 




~380 



Note: The figure of merit is computed as the non-normalized 
l/{(j{wp) X a{wa)) as in the DETF report. We have included the 
prior from Stage II dark energy experiments, prior on Hq from the 
Hubble space telescope and Planck prior. 

Table 2 

Comparison of figure of merits (FoM) 



BAO with BOSS is --71. Voids could thus yield very 
good constraints on Dark Energy from deep redshift sur- 
veys with smaller sky coverage. 

Finally, we note that in all this work we have only 
used the shape of the stacked voids as a probe. Addi- 
tionally, the re is inf ormation in the distribution of void 
Jones & Fry (1998). The distribution of the 



orientations 



intrinsic shapes of individual voids contains a great deal 
of complementary information, at t he potential cost of 



requiring to assume a galaxy bias ( Park fc Lee l 12007 
Biswas et al.||2010[ [Lavaux fc Wandelt||201U[ ).~We tHmk 
that the intrinsic shape could increase substantially the 
constraints that only comes from void shape analysis. 
This could be done in two steps: the first pass would use 
the method that we have developed in this paper to con- 
strain the local geometry, the second pass would use the 
shape distribution to constrain the growth of structures. 
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6. DISCUSSION AND CONCLUSION 

We showed that by identifying and stacking cosmic 
voids in redshift shells and size bins, and then measur- 
ing their shapes in redshift space we can directly con- 
strain the cosmological expansion through a purely ge- 
ometric approach. Several steps are required to use the 
stacked voids technique to connect a spectroscopic sur- 
vey to the expansion geometry and hence to dark energy 
phenomenology. In this paper we proposed methods for 
each one of these steps which, together, amount to a first 
analysis pipeline for the stack ed voids techniq ue. 

We use a modified ZoBOV ( Neyrinck]|2008 ) algorithm 
for finding and stacking voids on a light-cone, extended 
to produce non-overlapping voids selected according to 
two criteria: an effective radius within a given range and 
a central density sufficiently low to mark the region as a 
void. 

In Section [4j we applied this algorithm to mock light- 
cone catalogs obtained from three A'^-body simulations. 
We have tested the method in the original comoving co- 
ordinates of the simulation, which has provided us with a 
model of the density profile of the stacked void. Then, we 
have simulated cosmological expansion and distortions 
due to peculiar velocities, which allowed us to qualita- 
tively estimate the impact on our measurement of the 
stretching of voids. We find that even a crude de-biasing 
prescription of the (mild) peculiar velocity systematics 
yields a powerful method. By generating void stacks for 
different void sizes and redshift shells, we project out the 
details of individual void shapes. 

We have developed and tested a Gaussian statistical 
model able to estimate the stretching of the stacked 
voids. Combining the results obtained from void stacks 
at all redshift shells results in an estimate of the expan- 
sion history. 

We are aware that each one of these analysis steps can 
likely be improved significantly. There are parameters to 
optimize, such as the widths of the redshift shells, and 
the size bins. We only scratched the surface of possible 
methods in spatial statistics and computational geome- 
try when it comes to defining voids in realistic surveys. 
Our method for measuring the redshift space shapes of 
stacked voids is only one of many that one could imag- 
ine. In particular, we have not yet taken advantage of the 
ability to model the systematics due to peculiar veloci- 
ties, present at 10-15% of the expansion signal, except to 
suggest a crude de-biasing approach which gives a con- 
sistent result on our simulations. 

Based on these first results we extrapolated and per- 
formed Fishcr-matrix forecast of the constraints on Dark 
Energy equation of state we expect from the SDSS, BOSS 
and EUCLID spectroscopic surveys. 

We have found that cosmic voids have the potential to 
provide a far more powerful constraint on dark energy 
than measurements of the Baryonic Acoustic Oscillation 
scale, by an order of magnitude. This large increase of 
information is easily understood in comparing the num- 
ber of modes probed by voids compared to BAOs, which 
scales roughly as the third power of the ratio of the BAO 
scale to the scale of the smallest usable voids ~ 1000. 
The area of parameter constraints scales as the square 
root of the number of modes ~ 30. When projected into 
the Wa , Wp plane using the Fisher matrix formalism for 



the EUCLID wide survey, we find the improvement over 
BAO on those parameters by a factor of ~ 30. 

We expect our stacked void shape measurements to be 
robust to galaxy bias as it is purely geo metrical and re- 
lies o n the topology of the density field ( Springel et al. 
19981. In fact, it is possible that biased tracers ot the 
density enhance the contrast of voids and therefore en- 
hance the void detection rate. Another limitation comes 
from the effective volume accessible from galaxy surveys 
which is not infinite and may only cover some parts of 
the sky. If the survey area is sufficiently contiguous with 
respect to the considered void sizes then the void identifi- 
cation should not be affected by the constraints imposed 
by the geometry. Prior to the stacking, the identified 
voids must be rotated such that the stretching direction 
is always in the same fiducial direction. Finally, we do 
not expect apparent magnitude limitation to impair the 
analysis differently than just removing the smallest ob- 
servable voids at a given redshift. This expectation is 
based on the conservation of topological properties un- 
der resampling of the density field. All these expectations 
remains to verified on more realistic mock catalogs and 
real data. 

Based on our Fisher matrix forecasts, the stacked voids 
technique alone promises to double the figure of merit 
from EUCLID when compared to the combined results 
from all other probes using the same data (BAO, weak 
lensing, type la supernovae, cluster counts). The Alcock- 
Paczyhski test using stacked voids is therefore potentially 
a significant addition to the portfolio of major dark en- 
ergy probes which merits further detailed studies focused 
on additional real-world systematics and optimal survey 
design. 
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